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1 Introduction 


Recent years were quite productive in field of the particle spectroscopy. Many charmonium 
and bottomonium like states are observed by BaBar, Belle, LHCb and BES III collabora¬ 
tions [1]. These progresses in experiments stimulated further theoretical studies and ex¬ 
perimental investigations on this subject [2], Considerable progress has also been made on 
spectroscopy of the conventional heavy mesons states containing single charm and bottom 
quarks such as P,,(2700), T>:^(2860), T>,^(3040), T>j(2580), T>j(2740), T>}(2760), T>j(3000), 
P}(3000), i3i(5721), 132(5747), 8^,(5830), 13*^(5840), T>(5970), etc [3], Soon after DO Col¬ 
laboration observed the i3i(5721) and 132(5747) states [4], which were both conhrmed by 
the CDF Collaboration [5]. The CDF Collaboration further observed the Dsj(5830) and 
13*^(5840) states [6] which in turn conhrmed by the DO Collaboration [7]. Moreover, the 
masses of the (5830) and B*^ (5840) states were determined more accurately by the LHCb 
Collaboration [8]. 

The masses and decay constants of the heavy tensor meson D2(2740) and V*^{2573) 
states were hrst studied within the framework of the QCD sum rules method in [9,10]. 
These mesons, as well as the 132(5747), 13*^(5830) and 13*^(5840) tensor mesons have recently 
been studied within the same approach in [10]. Note that light tensor mesons without, and 
with strange quark have also been analyzed in framework of the QCD sum rules method in 
[11] and [12], respectively. 

One of the most promising approaches in investigating the properties of mesons and 
hadrons is the study of the electromagnetic form factors and multipole moments. These 
studies can provide useful information about their internal structures. The electromagnetic 
properties of usual mesons, as well as photons and neutrons have comprehensively been 
studied from theoretical and experimental sides, and at present it is the subject of the 
growing interest from both sides. However the study of the electromagnetic properties of 
the tensor mesons has received less interest, and therefore more effort is needed in this 
respect. The magnetic moments of the light tensor mesons were investigated in framework 
of the light cone QCD sum rules method (LCSR) in [13]. In the present work we calculate 
the magnetic dipole moments of the recently discovered heavy tensor mesons in LCSR. 

The paper is organized as follows. In section 2, the light cone QCD sum rules are 
constructed for the magnetic dipole moments of the heavy tensor mesons. In section 3, the 
numerical analysis is performed for the obtained sum rules. 


2 Theoretical framework 


Before presenting the light cone sum rules for the magnetic dipole moments of the heavy 
tensor mesons, let us hrst introduce the matrix element which corresponds to the transition 
of the heavy tensor meson with momentum p-|-g to the heavy tensor meson with momentum 
p in presence of the electromagnetic held, i.e.. 


{TQ{p,e) \jf\TQ{p + q,e))=el^{p)\2p^ 


g^Y^F,{Q^)-g 


ar 


2m^Q 




2 m^^ 




+ [rq'^ - g^q^) 




Fsm 

\£xT{p + q) ,( 1 ) 


1 



where Fi{Q‘^) are the form factors, and Tq{p, e) means heavy tensor meson with momentnm 
p and polarization tensor Eap- Since we are interested with the magnetic dipole moments of 
heavy tensor mesons, the valnes of the form factors at = —q^ = 0 need to be calculated. 
The transition under consideration can be described by the following correlation function, 

q)^-Jd‘xI (o \T{uM]t(y)Mx) }Io) , (2) 

where is the interpolating current of the heavy tensor meson, and is the electromag¬ 
netic current given as, 

fp = WpQ + egQ^pQ , 

with the electric charges and eg of the light and heavy quarks, respectively. The coupling 
of the tensor meson current to the tensor state is defined as, 

(0 \jpu\ Tq{p, e)) = ml^QTQefxu , (3) 

where rriTg is the mass, and qtq is the coupling constant of the tensor meson. 

It is more convenient to rewrite the correlator (2) by introducing the electromagnetic 
background field strength tensor 

of the plane wave, in the following form 

^pupai3£^ = i J (o T|j^^(a;)j„/3(0) 0^^ , (4) 

where F means that all vacuum expectation values are calculated in the background electro¬ 
magnetic held. The correlation function (2) can be obtained by expanding the correlation 
function (4) in powers of the background held, and taking into account the terms that are 
linear in F^^, which corresponds to the single photon radiation. The advantage of using the 
background background held is that it separates the hard and soft photon radiation in an 
explicitly gauge invariant way (for more about the details of the background held method, 
see [14] and [15]). 

After these preliminary remarks, we can now proceed deriving the light cone QCD sum 
rules for the magnetic dipole moments of the heavy tensor mesons. These sum rules can 
be obtained by calculating the correlator function in terms of mesons (physical part) from 
one side; and calculating the same correlation function in terms of quark-gluon degrees of 
freedom by using the operator product expansion (OPE) in deep Eucledian region from 
theoretical side. Matching these two representations of the same correlation function, the 
sum rules for the magnetic dipole moments of the heavy tensor mesons are obtained. 

Calculation of the correlation function from the physical side is performed by inserting 
the complete set of tensor meson states having the same quantum number as that of the 
interpolating current, and isolating he ground state, as the result of which we get. 


^pi/pafi 


(Tg(p, e) b'?I Tg(p + g, e)> ^ 
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where dots mean the contribution of the higher states and continuum. 

From experimental point of view the multipole form factors are more useful than the 
form factors given in Eq. (1). The relations between the two sets of the form factors for 
any arbitrary are derived in [16]. For the real photon = 0) these relations are given 
as: 


Fi(0) = G£;p(0) , 

-^ 2 ( 0 ) = Gmi(0) , 

-^ 3 ( 0 ) = —2GEoiO) + GE2i0) + Gmi(0) , 

^4(0) = —Ga4j(0) + Gm3(0) , 

-^ 5 ( 0 ) = G£;q( 0 ) — [G£; 2 ( 0 ) + G^Ari(O)] + G£;4(0) + 1^M3(0) , (6) 

where Ge^(O) and Gmi{0) are the electric and magnetic multipoles. 

Calculation of the correlator function from the physical side is performed by inserting 
the complete set of tensor meson states in Eqs. (1-4), from which for the hadronic part we 
get, 






{p‘^ - + qy - 


£f^Apyiaip)\ 2(p-£ 
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2m% 


2m% 2m% 






2m‘^Q 
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£xr{p + q)£a^{p + q) ■ (7) 


Summation over the polarizations of the heavy tensor mesons can be performed by using 
the relation, 


£*paiP)£Xr{p) = ^Fpxip)V^rip) + ^Vpr{p)V^x{p) - ^Vpaip)Vxrip) , ( 8 ) 

where 


Fpxip) = ( -9px + ^) 
^TqJ 


Using Eqs. (6) and (8) in Eq. (7), for the physical part of the correlation function we get, 
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It follows from Eq. (9) that the correlation function H^upafS contains various inde¬ 
pendent structures, each of which can be used in determination of the multipole moments 
of heavy tensor mesons. In the present work we restrict ourselves to the calculation of the 
magnetic dipole moment Gmi, and for this goal we choose the structure — e^q^)g^a- 
This structure has advantage over the others since it does not contain contributions coming 
from the contact terms (for more detail about the contact terms, see for example [17]). 
As the result we get the following sum rules for the magnetic dipole moment of the heavy 
tensor mesons. 


^pupajS ^ 9p( 






- ml )[{p + qf - ml 


Gmi + other structures 


.( 10 ) 


Denoting the coefficient of the {e^q’^ — e'^q^)g^a structure as the invariant function 11, we 
get the following hnal result for the physical part of the correlator function for the magnetic 
dipole moment Gmi, 


2 

‘^Tq9tq 1 

^ (p2 - ml^)[{p + qf - m%^] 

In order to construct the corresponding sum rules the correlation function needs 

to be calculated from the QCD side in terms of quark and gluon degrees of freedom using 
the operator product expansion, for which we need the interpolating current The 

interpolating current for the ground state heavy tensor meson with the quantum numbers 
2~^ can be chosen as. 



Jpiy 


q{x) 7 ^ Vv +'lu Vp Q{x) 


( 12 ) 


where q and Q denote the light and heavy quark helds, and the derivative operator Vp with 
respect to x^ can be written as. 



1 

2 


Vp {x)- Vp {x) 


which acts on the right and left sides, and the covariant derivatives in it are dehned as 


Vp {x) = dp {x) - f|AM“(a;) , 

Vp {x) = dp {x) + , (13) 

where A“ are the Gell-Mann matrices, and A“(a;) is the gluon held. In the present work 
we use the Fock-Schwinger gauge in which the external held A“(a;) satishes the condition 
a;^A“(a;) = 0. 

The correlation function is calculated from the QCD side in deep Eucledian region 
—)■ —oo, and (p -|- qY —)■ —oo after contracting the corresponding heavy and light quark 


4 



fields, as a result of which we get, 


(o \[s,{y - [ du (x) dy (y) 

- du {x) dy {y)- du {x) dy {y)+ du (x) dy iy) Sq{x - |/)7o| 

“h CX^ -\- {z/ i-X /ij" -j- {^13 i-X Q;, U i-X o\ , (1^) 

/ F 

where we set y = 0 after performing derivative with respect to y. 

We see from Eq. (14) that, the expressions of the propagators of the light and heavy 
quarks are needed in order to calculate the correlation function. The light quark propagator 
in presence of the external held is calculated in [18], whose expression is given as: 


s,(x - !,) = Sf“{x - 9) - M [l - i/)] + [l - ,/) 

- 'S, J ‘^“1 - y“)G„yu{x - y))Y 


where A is the scale parameter separating the perturbative and nonperturbative domains. 
This parameter is estimated in [19] to have the value A = (0.5-7 1.0) GeV] and S^''^'^{x — y) 
is the free quark operator whose expression is given as: 


S’yyx - y) 


Hf- i) 

2n‘^{x — yY 


mg 

Att‘^{x — yY 


( 16 ) 


The propagator for the heavy quark have the following form in coordinate space: 



where where Ki{mQ\/—x^) are the modihed Bessel functions, and 



Few words about the expression of the quark propagator are in order. The complete light 
cone expansion of the light quark propagator in presence of the external held is calculated 
in [18], which includes the contributions coming from nonlocal three qGq, and four-particle 
qqqq, qG^q operators. Using the expansion in conformal spin, one can show that aforemen¬ 
tioned contributions are small (for more detail see [20]), therefore we shall neglect them in 
further analysis. 
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There are three type of contributions to the correlation function: 1) Perturbative part, 
when photon interacts perturbatively with the quark propagator (light or heavy). 2) 
“Mixed” contributions, which take place when heavy quark propagator interacts with the 
photon held perturbatively, and light quark helds form quark condensate. 3) “Long distance 
” contribution. It takes place when photon is radiated at long distance. 

The perturbative contribution is calculated from Eq. (14) by replacing heavy or light 
propagator with, 

-v)^ , (18) 

where Tp are the full set of Dirac matrices. As has already been noted, when a photon 
interacts with the light quark helds matrix elements of the nonlocal operators such as 
q{x)Vq{ii) and q{x)Gqiy) appear between vacuum and photon states. These matrix 
elements are parametrized in terms of the photon distribution amplitudes (DAs), which are 
the key nonperturbative parameters in light cone sum rules, whose explicit expressions are 
given below. 


(7(g)|g(a;)ap^g(0)|0) = -ieq{qq){e^qy - e^q^) j 


2{qx 


:eg(gg) 


, ex\ ( ex\ 

qy — 

' qx J \ qx J 




(7(g)|g(a;)7pg(0)|0) = (sf, - q^,—'] [ {u) 

\ / Jo 

(7(g)|g(a;)7p75g(0)|0) = —egfs^ef.yape'^q^x^ 
{x{q)\q{x)gsGf,y{vx)q{0)\0) =-ieq{qq) {e^qy - Syq^) J 
{'y{q)\q{x)gsG^yi'y5{vx)q{0)\0) = -ieq{qq) {e^^qy - Syg^,) j 

{l{q)\q{x)gsG^y{vx)^o,im{^)\^) = eqf3'fqa{£fxq,x - Si^qij) J 

{l{q)\q{x)gsGf,y{vx)i'ycq{0)\0) = eqfz^qo,{£^,qy - £yq^q) j 
{l{q)\q{x)aafigsG^y{vx)q{Q)\Q) = eg{qq) | (y9au - 

- (^9fiu - 

- [su - q^-^ + q^l^S^ g/3 


qj 
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ex 


qx 

ex 


- U/3-9/3- - iq^^Xa +qaXf,) ] q^ 


qx 


qx 


“ 1 “ ^q qq ) Qua ijlvXa “ 1 “ QaXu') 1 qfi 


qx 




+ —- quX^j,){eo,q^ - efjqa) j 

+ —{qaXfi - qfiXa){e^qu - e^q^,) f \ , 

qx J 


(19) 


where is the leading twist-2, 'ip^{u), 'ip°'{u), A and V are the twist-3, and h^{u), A, 

Ti {i = I, 2, 3, 4) are the twist-4 photon DAs, and x is the magnetic susceptibility. The 
photon DAs are calculated in [15] and their explicit expressions are given in Appendix A. 
The measure Da,- is dehned as 


Da,- = / dan 


dan 


dag6{l — aq — aq — ag) 


Separating out the coefficient of the structure {e^q'^ — e'^q^)gg,a from the QCD and the 
phenomenological parts of the correlation function and equating them, we get the magnetic 
moments of the heavy tensor mesons. In order to suppress the contributions of the higher 
states and continuum, we apply double Borel transformation with respect to the variables 
and {p + q)"^- As the result of these calculations we obtain the following sum rules for 
the magnetic dipole moment of the heavy tensor mesons. 




eu{glG^)M‘^ [-2mlX2 + m^Xs) 


11527r2 L 


3456mfe7r^ 


{eu{g‘lG‘^) + 96ebmb7r‘^(uu)) 


+ 4:eu7r 


18mb{uu) A(mo) -h 2ji{h^) + Aj 2 {hy) + {g,G ){uu)xX 7 i^o) + 36/3^mb^/’“(Mo) 


1 

;ebmlM'^ (X 2 - mils) + 


euM^ 


+ 8mb{uu)xXy{uo) 


327 ^ 2 V- / ■ gg 

fs'r f6^“(wo) - 

2ebl2 + ^ebmlls AebmlX^ + 2e^,m^X4 + 2{eb- eAjmlX^ 
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- qmlM^ 

327r2 ’’ 

-,mb 


uu) + CuigiG/ — 4{uu)A{uo) + 4(5 — 4u){uu)ji{h^) 


6912M2 

A0{uu)j2{h^) +mb(^- 8mb{uu)xX'riuo) + fz^{-8>ii{i)'') + 6^/>“(mo) - 4^/>^(mo) -F ^“'(mq)) 


+ 


3456M6 


eu{gsG )ml ^{uu) [A{uo) + 2ji{h^) + Aj2{h ^)) -h 2/3^mfe^/>“(Mo) 
eu{gsG^)ml{uu)A{uo) + 


3456mb7r2 


,^A{{glG'^) - l8ml)Ti‘^{uu)A{uQ) 
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(20) 


+ {glG"^) - 3 ml + 7 r^(^ 4(2 + u){uu)ji{h^) + 16{uu)j2{h^) + mb{-12mb{uu)xy^'y{uo) 
+ /37[-12ji(^/’’') - 2 ^/>“(mo) + (-2 + u){A^/J^{uo) - ^“'(mo))]} 
where 


Ml 

“ Ml + Ml ’ 




MlMl 

HITmI ' 


The functions in {n = 1,2), and ji(/(n)) are dehned as: 


Ji(/(«)) 
J2 (/(«)) 
T 

-^n 


duf{u) , 





Uo)f{u) 





where sq is the continuum threshold. 

Since we have the same heavy tensor mesons in the initial and hnal states, we can set 
Ml = Ml = 2M^, as the result of which we have. 


Ml 

“ {Ml + Ml) 


1 

2 ■ 


Physically this means that each quark and antiquark carries the half the photon momentum. 


3 Numerical analysis 

This section is devoted to the numerical analysis of the sum rules for the magnetic dipole 
moments of the heavy tensor mesons obtained in the previous section. The values of the 
input parameters entering into sum rules are, {uu){fj, = 1 GeV) = {dd){fi = 1 GeV) = 
— (0.243)^ GeV^, (ss)|^=i cev = {0.8±0.2){uu){jj, = 1 GeV), ml = (0.8 ±0.2) GeV^ which 
are obtained from the mass sum rule analysis for the light baryons [21,22], and B, B* 
mesons [23]. For the heavy quark masses we have used their MS values, which are given 
as: Thbirhb) = (4.16 ± 0.03) GeV and rhcirhc) = (1.28 ± 0.03) GeV [24,25]. The magnetic 
susceptibility of quarks was estimated in [26-28] in framework of the QCD sum rules. As 
we have noted earlier, the residues and masses of the heavy tensor mesons were calculated 
within the QCD sum rules method in [9,10], and their values are gTi 2 = 0.0228 ± 0.0068, 
gvs^ = 0.023 ± 0.011, gs^ = 0.0050 ± 0.0005, gss^ = 0.0060 ± 0.0005. 

Having decided the values of the input parameters, we are ready now to perform nu¬ 
merical analysis of the sum rules for the magnetic dipole moment of the heavy tensor 
mesons.The sum rule contains two unphysical parameters: a) The Borel mass parameter 
M^, and continuum threshold sq. It is known that the physical results should be indepen¬ 
dent of the these parameters. Therefore, our primary goal is to hnd such domain of these 
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parameters for which the magnetic dipole moment is practically independent of them. The 
“working region” of is determined as follows: The upper bound of can be found 
by requiring that the contributions coming from higher states constitutes about 40% of the 
perturbative part. The lower bound of could be hxed by demanding that the higher 
twist contributions are less than that of the leading twist contributions. In other words, 
the light cone expansion with increasing twist should be convergent. These requirements 
leads us to the following domains for the Borel mass parameter: 

2.0 GeV^ < < 4.0 Geld^, for V 2 and Vs 2 mesons , 

4.5 GeV^ < < 7.0 GeV^^ for B 2 and Vb 2 mesons . 

The second parameter entering to the sum rules is the continuum threshold sq. Generally 
speaking, this parameter is not arbitrary and it is related to the energy of the first excited 
state. The energy necessary for the transition of the meson from ground state to first excited 
state is yGo — niground- This difference varies, usually, from 0.3 GeV to 0.8 GeV, where in 
our numerical analysis we have used their average value, i.e., ^/so — niground = 0.5 GeV. 

In Figs. (1) ,(2) and (3) we present the dependence of the magnetic dipole moments of 
Pj, P 2 Vs 2 tensor mesons on M^, at various hxed values of the continuum threshold. 
From these hgures we get the following results, 

( 0.75 ± 0.25, for 

GmAq^ = 0) = < -2.10 ± 0.20, for 

[ -2.20 ± 0.20, for V~^ 


in units of e/2mTQ. 

Similar analysis performed for heavy tensor mesons containing 5-quark, whose results 
are presented in Figs. (4)-(6), respectively. From the analysis of these hgures we obtain, 

( 3.8 ± 0.7, for 

GmAq^ = 0) = < -1.3 ± 0.3, for B^ 

[ -1.4 ± 0.3, for B^S 2 


in units of e/2mTQ. 

In conclusion. The magnetic dipole moments of the heavy tensor mesons are calculated 
in framework of the light cone QCD sum rules. It is observed that the magnetic moments of 
the charged tensor mesons are larger compared to the neutral ones. The SU (3) symmetry 
breaking in heavy tensor mesons containing beauty quarks is about 10%, while in the 
charmed meson sector it is quite large. 
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Appendix A: Photon distribution amplitudes 

Explicit forms of the photon DAs [15]: 
r - 

(P'yiu) = 6uu 1 + (p 2 {fi)C 2 {u — it) , 

= 3[3(2m - 1)2 - 1 ] + ^(15< - 5w^)[3 - 30(2m - l)^ + 35(2 m - 1)^] , 

<(m) = [1 - (2m - 1)2] [5(2m - 1)2 - 1]^ 1^1 + ^ , 

A{ai) = 360Q;qQ;qQ;2 i + - 3) , 

V(aj) = 540w):^{ag - ag)agagag , 

h^{u) = — 10(1 + 2 k^)C2 {u — u) , 

A(m) = 40M2-u^(3fi; — + 1) + 8 (^ 2 ^ — 3 C 2 )['Wm (2 + 13m'u) + 2m^(10 — 15m + 6m 2) ln(M) 

+ 2m3(10 - 15m + 6m 2) ln(M)] , 

120(3((^2 a C 2 )(*^g Oig)OigOigOig , 

7i(ai) = 30a2(ag - ag)[{K - k+) + (Ci - C^)(l - 2 ag) + ( 2(3 - Aug)] , 

7i(ai) = - 120 ( 3(2 - C2^)(ag - Oig)aqagag , 

'T^ioii) = 30Q;2(Q;q — ag)[{K + K^) + ((1 + Cl")(l — 20 ;^) + ( 2(3 — AUg)] , 

~ 30Q;2{(ft; -|- fi;''')(l — ag) + ((1 + Ci~)(l — 0 'g){l — 2ag) 

A~ C2[3(cig ^q) Cl(j(l *^(;)]} ) 

§{ai) = -30a2{(M - «+)(! - ag) + (Ci - Ci+)(1 - ag){l - 2ag) 

A~ ( 2 ( 3 ( 0 :^ ^q) C^g(l *^<;)]}' 

The parameters entering the above DA’s are borrowed from [15] whose values are <^ 2(1 GeV) = 
0, w]^ = 3.8 ± 1.8, = -2.1 ± 1.0, K = 0.2, k+ = 0, (i = 0.4, (2 = 0.3, Ci+ = 0, and 

C2^ = 0. 
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Figure captions 

Fig. (1) The dependence of the magnetic dipole moment of the T >2 tensor meson on 
at two fixed values of sq = 9-0 GeV"^, and sq = 9.5 GeV^. 

Fig. (2) The same as Fig. (1), but for the T >2 tensor meson. 

Fig. (3) The same as Fig. (1), but for the tensor meson. 

Fig. (4) The same as Fig. (1), but for the tensor meson, at two hxed values of 
So = 37.5 GeV‘^, and Sq = 40.0 GeV^. 

Fig. (5) The same as Fig. (4), but for the tensor meson. 

Fig. (6) The same as Fig. (4), but for the 0^2 tensor meson. 
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